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The correct use of energy-dependent single-particle level (s.p.l.) densities 



(ESM) is analysed. First, an analytical expression is obtained following the 
convolution of energy-dependent excited-particle and hole densities. Next, a 
comparison is made with results of the ESM formula using average s.p.l. den- 

a ' 

sities for the excited particles and holes, respectively. The Fermi-gas model 
• <— i . 

(FGM) s.p.l. densities calculated at the corresponding average excitation 



energies are used in both cases. The analysis concerns also the density of 
particle- hole bound states. The pairing correlations are taken into account 
while the comparison of various effects includes the exact correction for the 
Pauli exclusion principle. Quantum-mechanical s.p.l. densities and the con- 
tinuum effect can also match a corresponding FGM formula, suitable for use 
within the average energy-dependent partial state density in multistep reac- 
tion models. 
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I. INTRODUCTION 



The particle-hole state densities are basic quantities for the description of preequilibrium 
emission (PE) in semiclassical models as well as quantum-statistical theories (e.g., [0,0) 
involving a series of particle-hole excitations caused by two-body interactions. The nuclear 
excitation in the equilibrium processes concerns the single-particle levels (s.p.l.) within an 
energy range of the order of the nuclear temperature around the Fermi level. This explains 
the basic role of the s.p.l. equidistant spacing model (ESM) || in the analysis of the 
equilibrium emission (see also M). However, much higher and lower single-particle energies 
are involved in PE reactions so that one should consider the reduced suitability of the ESM 
partial-state density (PSD) formula of Williams ||. Moreover, the inconsistency between 
the phenomenological s.p.l. density g ~ A/14 MeV" 1 and the number A of nucleons in 
the nucleus has come under increasing criticism ||[7|. On the other hand, combinatorial 
calculations performed in the space of realistic shell model s.p.l. || have other inherent 
shortcomings (e.g., the strong dependence on the basic set of s.p.l.) HJlO[|. This accounts 



for the continued use of the Williams- type formula with various corrections [11,12] or exact 



calculation [13[ for additional Pauli blocking and the pairing interaction. 

In fact, there have been early attempts at considering the single-particle energy de- 
pendence of the s.p.l. density g(s) within PE formalisms [H-|nj]. Next, Kalbach [ JTT|JT5[ 
discussed different forms of this dependence and found it tied to PE surface effects due to 
the interdependence of the respective assumptions. Herman et al. || obtained an indication 
for the energy dependence nearly as that given by the Fermi-gas model (FGM) below the 
Fermi energy F, but linear above F. Chadwick and Reffo |7j] found the use of either the 
FGM prescription or the equidistant parametrization g=A/F more accurate than the phe- 
nomenological one. The FGM s.p.l. density has also been involved in the development of 



the partial level densities with linear momentum (T£) . At the same time, the ESM accuracy 
has been discussed in connection with the non-uniform s.p.l. density effect |2(| provided by 
the harmonic oscillator model. The analysis of the energy-dependent s.p.l. density in the 
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vicinity of the Fermi energy [21j provided a more general form and a good approximation of 
the effect for low energies, where the influence of the finite depth of the potential well can 
be neglected. Various g(e) have been obtained within both the semiclassical Thomas- Fermi 



approximation [pl|-p8| and the exact quantum mechanical calculations p9|-pT| which are 
also applicable at the high excitations located in the continuum region. The PSD includ- 
ing distinct energy-dependences for the excited-particle and hole level densities has recently 
been used in semiclassical []32| or quantum-statistical f33|,|34| cross-section calculations. 

The valid use of energy- dependent s.p.l. densities within the ESM particle-hole state 
density formula, even when corrected for the finite depth of the real nuclear potential well 
||35|| , has not yet been proved. Proving it is one aim of this work. First, the particle-hole state 
density is obtained in Sec. II by means of recursive relations particularly using the FGM 
s.p.l. density. Next, these are compared in Sec. Ill with the results of the ESM formula 
modified by using s.p.l. densities different for excited particles and holes, obtained from 
the FGM at the respective average-excitation energies [[OJ (the average energy-dependent 
ESM formalism). The analysis is also carried out for the density of particle- hole bound 



states, with single-particle excitations not exceeding the nucleon binding energy [3H]. The 
advanced pairing correction JTT| , [r2| ] is taken into account while the comparison of various 
effects includes the exact correction for the Pauli exclusion principle [fl3 1 . The importance of 
distinct corrections in the average energy- dependent ESM formalism is further discussed in 
Sec. IV. At the same time the subtraction of the free-gas contribution [[2^,[3(|[J7| is analysed 
within this formalism, thus making no use of arbitrary truncation [|38[| . The respective results 
are compared with the semiclassical and quantum-mechanical calculations of the continuum 
effect. Since the actual quantum-statistical analyses of the multistep reactions still involve 
the rough ESM, the respective results could be altered following consideration of the effective 
iViV-interaction strength as the only free parameter. The conclusions are drawn in Sec. V. 



II. THE PSD RECURSIVE FORMULA 
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A. Single-particle level densities 



Densities of the excited particles and holes with distinct energy dependences or even 



different values at the Fermi energy F were considered by Gadioli and co-workers [15,16 
Betak and Dobes ||35|| , and Herman et al. ||39|| . The subsequent study of unperturbed 
shell-model Hamiltonian spacings indicated a linear energy dependence for excited particles, 
as well as different corresponding values at the Fermi level. On the other hand, Schmidt 
et al. EDI found that the smooth s.p.l. density in a Woods-Saxon potential lies between 



the density corresponding to an infinite box and the one for an harmonic oscillator, and 
approximately follows g(e) ~ e. Moreover, this energy dependence has already been used 
within an improved abrasion model for heavy-ion collisions JIT] . 



Given the need for an analytical PLD expression, we have followed the method of Bogila 



et al. |21| while the finite depth of the nuclear potential well and the case of particle-hole 
bound states have also been considered. Actually, the particle-hole bound state formula 
turns into the common form in the limit of large values of the nucleon binding energy B. 
The following discussion will concern the general form of g{e), with the go value at the Fermi 
level. However, the usual FGM energy dependence 

is particularly taken into account (see also the Appendix). This can be expressed in terms of 
the single-particle excitation energies u = e — F for particles, and u = F — e for holes. Next, 
similarly to Bogila et al. we have retained the first three terms of its expansion in powers 
of u around the value at zero excitation energy. The general forms of the excitation-energy 
dependence then become 

g p (u) = (au 2 + bu + c)6(B - u) (2) 

and 

g h (u) = [au 2 -bu + c)6(F - u) , (3) 



where the theta functions are unity if their argument is greater than zero, and zero otherwise. 
The FGM values of the coefficients are a = = —g /8F 2 , b = g' = g /2F, and c = g , 
where g' and g' ' are the values at the Fermi level of the respective derivatives. Various 
energy dependences of the s.p.l. density can obviously be involved within this framework, 
by using the appropriate values for the coefficients in Eqs. (0-0). 



B. The convolution state-density formula 

The bound-state density uj(p,h, E) for p excited particles above the Fermi level and h 
holes below it (n = p + h), at the total excitation energy E, can be obtained by convolution 
of the single-particle and hole level densities with an excitation-energy conserving delta 
function 0,00 

u(p,h,E) 



(4) 



d«i0p(iti) / du 2 g p {u 2 ) ... / du p g p (u p ) 
p\h\ Jo Jo Jo 



x / du x g h {ux) . . . 



/ du h g h (u h )5 \ E -J2u x -J2 u j) > 



where the Pauli principle is not yet taken into account. One way to proceed || is to replace 
the 5 function by its integral representation 
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so that 



w(p, h, E) 



e ikE 



g p {u)e iku du 



(6) 



2vrp!/i! 

By using the s.p.l. densities given by Eqs. (|2|-|3]), evaluation of the integrals, and expansion 
of the respective results, it results 

exp[ik(E - XB-jF)] 



u{p, h,E) 



9o 
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where, by replacing A by z, 



p-i ki h-j h / r/\h+h /„/\ fc i+ z i 
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(So)* +J 



x + fl&S + J 2 QtfS + ^) fe2 - i2 - g' F + g J ' «F - <^ 2 ^ z , (8a) 

z is the integral which forms the functional argument in Eq. (0), and 

N = n + h + jt + h + h + i 2 + 32 + k 2 + h ■ (8b) 
Finally, by using the Cauchy residue theorem we have obtained the expression 

n P h 

w fo E) = ihU -iM S Ei-^CpCUE ~ iB ~ 3F) n - X 0{E -iB- JF) 
p.n.\n i). i=0 j =0 

x^((^-^-J^-"^|) ■ 0) 
The (n — 1) factorial and power have additionally been included in order to obtain a form 



similar to those obtained by Oblozinsky |36[ in the frame of the ESM, and Bogila et al. [2T 
within the FGM with no constraints for particles or holes. Thus, in the limiting case of large 
B and F in the above expression, the % and j indices become zero so that the last four sums 
in Rij disappear and the formula of Bogila et al. is obtained 

u( P ,h,E) = u; E (p,h,E) E E E EC-i^cfc&cfrcf 

fc 1= o h=aii=oji=o 



\9o) \9o) ' (n-l + h+jx + h + h)! ' 

where 

p\n\[n — 1)! 



x 
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is the well-known Ericson formula || for the ESM case. However, we underline that a definite 
single-particle ground state should be marked out by the finite value of F within a consistent 
energy-dependent s.p.l. density. Therefore, the finite depth of the nuclear-potential well 
should be explicitly present in the particle-hole state density formulas exceeding the ESM 
framework. The usual formula for the density of the particle-hole bound states within the 
ESM approximation results immediately from Eq. @ by noting the unity values of all 
functionals Rij for a constant s.p.l. density. The only difference with respect to the formula 
of Oblozinsky [|36| concerns the Pauli-blocking factor A p h and the corresponding minimum 



energy a p h for a p-h state, which have not yet been included here. On the other hand, the 
temporary omission of the Pauli-principle correction can be used to estimate the energy- 
dependence effect better. Thus, the ratio between the results of Eq. (|9D and Oblozinsky 
formula without the Pauli correction is shown in Fig. 1 for a few simple p-h configurations. 
Here the bound-state condition is released, <7o=14 MeV -1 , and F=38 MeV. These particle- 
hole states are actually the most important ones for the PE description, and are least affected 
by the Pauli-principle oversight. Also shown is the ratio of the PSDs given by Eq. (|9[) and 
the Ericson formula, in agreement with the trend of the Bogila et al. results (Fig. 1 of Ref. 
|21|| ) if one takes into account the different F values used in these analyses. While the former 
ratio describes the energy- dependence effect versus an ESM formula including the potential 
finite-depth, the latter has the same role versus the simplest ESM expression of Ericson. 

First, the case of the lplh configuration, Fig. 1(a), shows that the deviation of gh{u) 
from g is not really compensated by the corresponding deviation of g p (u), within the s.p.l. 
density convolution. Thus, the Oblozinsky and Ericson formulas provide the same PSD 
at excitation energies below F, which are higher than the results of Eq. ([|). Second, a 
consideration of the potential finite depth decreases the Oblozinsky PSD values versus the 
Ericson formula. At the same time, the results of Eq. (|9|) decrease less significantly above 
E=F because of the deviation of g p (u) from go. Moreover, the case of the 2plh configuration 
shows the increased importance of the energy dependence versus the Oblozinsky formula, 
with the deviation from unity of the respective ratio becoming relevant at excitations higher 



than F. Smaller deviations have been obtained 
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as compared with the Ericson expression. 



C. Inclusion of Pauli blocking and pairing correction 

Following the related forms of the recursive PSD expression (|9|) and Oblozinsky formula 
the correction for the Pauli blocking and pairing effects can be implemented within 
the former by inclusion of (i) a p-h configuration-dependent threshold energy in the theta 
function, and (ii) the Pauli-blocking and pairing-correction term of the excitation energy, 



within the Kalbach formulation [|T^] 



p(p+l) + h(h + l) (p - l) 2 + (h - l) 2 
Ak(p, h) = E th (p, h) ; + 



4#o 



g F(p, h) 



(12) 



where the threshold energy 



E th\P, h ) = a + Pm 



Pm\ +A 2 

,9o 



1/2 



(13) 



is determined by the ground and excited-state gaps A and A(p,h,E). The ground-state 
gap is related to the condensation energy C = g^A^/A which can be given by the constant- 
pairing correction U p [pUf , based on the odd-even mass differences (e.g., 0]). A is obtained 
by using the parametrization [|TT|,0 



A = [0.996 - 1.76(n/n c ) im (E/C)- 0S8 } 6(E - E phase ) , 



(14) 



where n c = 0.792g Ao is the critical number of excitons, and E p h ase is the energy of the 
pairing phase transition given by 



E phase = C [0.716 + 2M(n/n t 



\ 2.171 



(n/n c - 0.446) . 



(15) 



Actually, the latter theta function has been introduced by Kalbach [0 in order to explicitly 
take into account the lack of a phase transition for small n. 

The inclusion of p m =maximum(p, h) and the form of the second term in the Kalbach 



correction (|TJ) have been adopted for a Pauli-correction function symmetric in particles and 



S 



holes, including the effects of passive holes. Next, the third term in Eq. (^2|) has been 
added in order to force the PSD to have the values of g and 2g for E = E^(p, h) and 
E = Efaip, h) + l/<?o> respectively. The function 

F(p, h) = 12 + Ag [E - E th (p, h)]/p m (16) 

restricts the action of this third term to just around E^(p, h). 

Consequently, the PSD recursive formula (|J), now including the Pauli and pairing cor- 
rections, becomes 

p h 



U <P> K E) = n iftir„-lV S U-^ +J C;C{[E - A K (p, h)-iB- jF]- 1 

x 9{E - E th -iB- jF) J2y ([E - A K (p, h) - iB - J^^^j ■ (17) 
The completeness of the Eqs. (Q) or ( |T7D has the drawback of making them difficult to 



use in reaction calculations, due to the intricate form ( 8aj) of the functionals R 



III. THE AVERAGE ENERGY-DEPENDENT FORMULA 

Since approximate but simpler solutions are still of real interest, we will discuss below the 
Kalbach fl8j attempt to use the energy-dependent s.p.l. density within the ESM formula. In 
fact, we will be checking its correctness against the exact expression flT7|). It seems important 
to note that the Kalbach approach involves distinct but average s.p.l. densities for the holes 
and excited particles, respectively, at average excitation energies. 

A. The finite-depth and pairing corrections 

The general form of the ESM density of particle-hole bound states corrected for (i) the 
Pauli exclusion principle ||, (ii) pairing interactions 0,0, and (iii) the finite depth of the 
nuclear potential well p5||36| can be written, similarly to Kalbach fll8| , 

«(p, h, E) = —f K (p, h, E, F) , (18) 

p\n\(n — 1)! 
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where 



i=0 j=0 



E-A K (p,h)-tB-jF 1, - > 



E 

x9(E-E th -iB-jF) (19) 

is a function including the finite-depth, Pauli-blocking and pairing corrections, as well as the 
bound-state condition. The original correction function [pl| was used to modify the PSD 
formula for the infinite potential well, with the Pauli correction terms A(p, h) neglected in all 
terms except the leading one. However, by means of Eq. ( p~8|) the function fx(p, h,E,F) can 
now be regarded as the ratio between the actual PSD formula and the Ericson expression. 



The index K is related to the inclusion of the advanced pairing correction of Kalbach [|12 
i.e., of the terms Ak(p, h) and -E^lr 

The importance of the Pauli-blocking and pairing correction term A^(p, h) is shown in 
Fig. 2 for simple p-h configurations. First, omission of this correction provides the unit 
value of the function fx in the case of the xplh configurations at excitation energies lower 
than F. Second, there is a rather similar case (or even identical for the state lplh) if 



the Pauli correction A v \ t and the corresponding minimum energy a p h |36| are taken into 
account but not the pairing corrections. A small threshold behavior becomes apparent in 
this case. Third, the inclusion of the advanced pairing correction yields a strong reduction at 
lower excitation energies. Finally, the bound-state condition obviously provides a different 
function, mainly determined by the number of holes. 



B. Single-particle average excitation energies 

In order to take into account the long-range deviations from ESM, Kalbach (18| proposed 
the use of average values for the FGM s.p.l. density (fj) corresponding to average excitation 
energies for either particles or holes. As a first approximation, these energies were estimated 
in the ESM frame. We follow the same method below but also include the case of the bound 
states. 
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The probability for the occurrence of ap-h state with the excitation energy E and an ex- 
cited particle between u and u+du is given by g p (F+u) -uj(p — l, h, E — u, F)du/co(p, h, E, F) 
Consequently, the average excitation energy per excited particle is given by 

1 r B u- g p (F + u)u(p-l,h,E -u,F) 



u. 



du 



(20) 



p Jo u>(p,h,E,F) 

where B = minimum(E , B) . Assuming a slow energy dependence of the correction term 
Ak, the average excitation energies for either particles or holes become 

Ef+(p,h,E,F) 



u. 



n Jk(p, h,E,F) 



u h 



E — pu p 
h 



(21a) 



(21b) 



where 



f+(p,h,E,F) = tE(-l) ,+3 ^ 

i=0 j=0 



E-A K {p, K)-iB-jF 



E 



x 



1 + 



?? 



%B 



P E-A K {p, h)-iB-jF 
returns to fx(p+ 1, h, E, F) for large B 



9(E-E th -iB-jF) 



(22) 



The shapes of the two functions fx(p, h, E, F) 
and fx(p,h,E,F), and the average excitation energies for simple p-h configurations are 
shown, bound-state case included, in Fig. 3. The above-mentioned similarity between the 
functions fx(p, h, E, F) and fxip + 1, h, E, F) can be observed, in the general case, for the 
Iplh and 2p2/i-states in Figs. 3(a) and 3(b). 

It is worth noting the results of Eqs. ( |21| ) for the bound-plus-continuum states shown in 
Fig. 3(e), and the bound states only as displayed in Fig. 3(f). The distinct trends are due 
to the separate constraints on hole excitation up to a value of F in the former circumstance, 
and a particle excitation limited by the B value in the latter. First, u p and Uh increase nearly 
as E/n at the lowest excitations, the slightly larger values for holes arising from the fact 
that fx(p,h,E,F) is smaller than fxip,h, E, F). Next, the constrained average excitation 
energy of either holes (in the general case) or particles (for bound states) becomes rather 
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saturated at energies E above the values of F and B, respectively. The ESM basis of Eqs. 



(pl|) determines the saturation values around F/2 for u h in the former case, and around B/2 
for u p in the latter. Moreover, in the limit of the lplh configuration a quite sudden change 
can be observed at total excitation energies around the values of F and B, respectively. 
There is also a small change in the trend of Uh and u p just below the maximum excitation 
energy E for a given p-h bound-state configuration, due to the final occupation of the highest 
allowed single-particle levels. 

To underscore the correlation between the specific shapes of the correction functions 
fx, average excitation energies u p and Uh, and corresponding average values of the s.p.l. 
densities 

9 P (p,h) = g(F + u p ) (23a) 

g h {p,h) = g{F-u h ), (23b) 

these quantities are shown together in Fig. 4 for the configuration 2plh. First, the effect of 
the Pauli and pairing correction is rather small. The saturations of u p and Uh are distinctly 
caused by the bound-state condition, Fig. 4(c), and finite-depth correction, Fig. 4(d), 
respectively. This is why the energy dependence of g p (p,h), Fig. 4(e), is distinct from the 
one of gh(p,h), Fig. 4(f), either in the general case (solid curves) or for the particle-hole 
bound states (dotted-dashed curves). 

C. ESM formula with average energy-dependent s.p.l. densities 

The average energy-dependent formula was finally obtained by using the average s.p.l. 
densities ( |23"D within the PSD formula ([T^) which becomes 

^Q= ^^W*- 1 f K{p ,H,B, F) . (24) 

It approximately takes into account the energy dependence of the s.p.l. density, even though 
the simple ESM form is still in use. However, there is no basic argument why this procedure 
should be used so that its accuracy needs further study. 

12 



The method adopted in this respect consists in a comparison of the results obtained 
by means of the approximate formula and the recursive Eq. fllTf). The corresponding 
predictions and their ratio are shown in Fig. 5 for the same simple p-h configurations. The 
global values F=38 MeV, £=10 MeV, g =U MeV -1 , and A =l MeV were used. First of all, 
in the general case of the particle-hole state densities (i.e., for B — > oo) there is only a small 
difference between the two PSD formulas even at medium energies. The agreement improves 
for more complex configurations, where the average of the single-particle excitation energies 
becomes really meaningful. Similar agreement is seen when the respective densities for the 
particle-hole bound states are compared within the first half of the energy range for each p-h 
configuration. However, the difference becomes significant near the maximum energy for a 
given particle-hole state, and increases as the PSD values drop back to zero. Nevertheless, 
the disagreement of these bound-state density formulas at the high-energy extremity should 
have little or no effect on the reaction cross-section calculations. 

Therefore we may conclude that the results obtained by using the average energy- 
dependent s.p.l. densities within the ESM formula are rather close to the exact convolution 
of the energy-dependent s.p.l. density. The next question concerns the need for this av- 
erage energy- dependent approach. The answer can be obtained by comparing the average 



energy-dependent ESM results with the PSDs given by the widely used ESM formula |36 
This is shown in Fig. 6, where the above global parameters were used. The pairing effect 
is apparent within this latter comparison, especially at lower energies. Besides this aspect, 
the behavior shown at medium energy is similar to the comparison of the recursive and 
Oblozinsky formulas in Fig. 1. The overall difference obviously exceeds the variation be- 
tween the predictions of the PSD recursive formula and the average energy- dependent ESM 
formalism. It is rather small for the general case of the particle- hole states, Figs. 6(a) and 
6(c), but larger for the bound states, Figs. 6(b) and 6(d). This strong effect following the 
consideration of the energy- dependent s.p.l. density is particularly caused by the constant 
increase in hole excitation for larger E, and the related significant decrease of the hole-state 
density as shown in Fig. 4(f). 

13 



D. Effect of exact Pauli-correction calculation 



The Pauli-exclusion effect on particle-hole state densities has already been subject to 
additional investigation by Zhang and Yang who used an exact method. Kalbach estab- 



lished later [12] that no conflict exists between their results and the frequently used Williams 
formulas if the energy-dependent Pauli term included pairing and passive-hole effects. The 
ESM derivation of PSD formulas without any approximation in the Pauli correction term 
was performed by Baguer et al. and Mao Ming De and Guo Hua The latter ex- 
tended the method to the case of the finite-depth potential and bound states, and included 



the Kalbach |L2[ pairing correction. The effect of the alternative use of the approximate and 
exact Pauli-correction, and the one caused by the average energy-dependent ESM formula 
are compared below. 

The results that were obtained by means of Eq. ( |18D and according to the exact Pauli- 
correction formalism |13[ are shown together in Fig. 7, for the most sensitive low-energy 
region. The analysis was first carried out without pairing correction, i.e., for Ao=0 as 
shown in Figs. 7(a) and 7(c). The global value go=14 MeV -1 adopted by Fu and Kalbach 



11,12] was used. Then, the PSD calculation with pairing correction corresponding to the 



value A =l MeV, Figs. 7(b) and 7(d), completed the analysis of these effects under any 
circumstances. 

The ratio of the PSD obtained by using the two formalisms, Figs. 7(c) and 7(d), shows 
more exactly that a close agreement - even within 1 % - is established just above the thresh- 
old for each p-h configuration. Minor deviations exist only for larger number of excitons, of 
less interest for multistep reaction calculations which include them in the so-called r-stage 
@. Therefore, the results of the exact Pauli-blocking effect calculations are closely related 
to those obtained by using the approximate Pauli correction |5|,[^,[36|] . The inclusion of a 
suitable pairing correction seems more significant with the following additional remark. The 
PSD for the very-few-exciton configurations become rather saturated within a short energy 
range above the threshold. On the other hand, these configurations have the main role 
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in the description of multistep reactions. Thus, the adequate account of the pairing effect 
may be found unessential for some analyses. The analysis of the high-energy limit of the 
particle-emission spectra, however, is quite sensitive to both pairing and nuclear-shell effects 
(e.g., ©0,0)- 

Finally, the comparative analysis of the effects illustrated in Figs. 6 and 7 shows, on 
a common basis, the higher importance of the s.p.l.-density energy dependence versus the 
exact calculation of the Pauli correction. 



IV. REALISTIC AND GLOBAL RESULTS 

The importance of the various approximations involved in the derivation of the PSD 
formulas should be well known in order to avoid useless effort or deficient results. First, one 
might want to know the consequences of the present analysis on the total state densities 
obtained as the sum of all PSD for allowed particle-hole numbers p—h. Next, one might 
question the usefulness of the PSD formulas discussed above, while quantum-mechanical 
calculations concentrating on the continuum region are being developed pipi . 



A. Effect on total state density 

Average s.p. excitation energies are shown in Fig. 8 for representative p-h configurations 
of the PSDs sum defining the total state density u(E). The same global values go=14 
MeV -1 , F=38 MeV and A =l MeV were used as above. Thus, it becomes apparent that all 
significant terms of this sum are characterized by rather equal u p and Uh increasing nearly as 
E/n. The saturated average excitation energy of the holes plays a major role in the case of 
the few-exciton states which are vital for PE description, but not for configurations around 
the most probable exciton number n |[46|| . The latter class of configurations mainly determine 
the total state density value, so that the corresponding ESM predictions are meaningful in 
this respect. Nevertheless, the average energy-dependent approach should be considered for 
the calculation of the PSD involved in the first stages of the multistep processes. 
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B. Consideration of the continuum effect 



Shlomo |3(J performed an exact quantum mechanical calculation of the s.p.l. density 
as the sum of the bound and continuum contributions in the case of finite potential wells. 
A distinct point of this approach has been the consideration of the free-gas states counted 
by the s.p.l. density for a finite potential well. The density of these states was calculated 
and subtracted by using Green's functions associated with the respective single-particle 
Hamiltonians. Then, the commonly used semiclassical approximations for the s.p.l. density 
were similarly considered for some widely used mean-field potentials. Thus, Shlomo found 
by means of both classes of methods that, for a realistic finite depth well, the s.p.l. density 
decreases with energy in the continuum region (the continuum effect). 

This effect may have a twofold meaning for the multistep reaction calculations. First, 
the continuum s.p.l. density following the subtraction of the free-gas contribution should be 
added to the particle-hole bound state density. The latter quantity has been used for the 



description of the multistep compound (MSC) processes J47j. It has been assumed to be 



zero outside the nuclear well, which is now considered less appropriate f2~7 |. Second, in the 



opinion of Bogila et al. |28| the subtraction of the free gas spectrum should be involved in all 
PSD for PE calculations. This point could be most important in accounting for the multistep 
direct (MSD) processes [37| which currently take into account all particle-hole states. Thus 



a correct yet simple method to estimate the s.p.l. density including the continuum effect is 
needed. 

The continuum effect can be taken into account within the average energy- dependent 
ESM formula (|4]) by using a form similar to Eq. (30) of Ref. (see the Appendix) for 
the excited-particle level density. According to Eq. (|23a|) it becomes 



9 P (P, h) = |^ 



F 



1/2 /<7-mV2 



l + I ~~r- J 6(u v -B 



"p 



(25) 



To emphasize the origin of the particular behaviour of g p (p,h), this is shown at the same 
time as the corresponding u p for the basic 2plh configuration in Figs. 9(b) and 9(d). It is 
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obvious that the average excitation energy is unchanged (see Fig. 4) so that the continuum 
effect fully determines the corrected s.p.l. density for excited particles. The comparison with 
the similar quantities Uh and gu(p, h) (the same as in Fig. 4) demonstrates the importance 
of this effect, the average s.p.l. density becoming even lower for excited particles than for 
holes. 

On the other hand, it seems worth comparing the specific average excitation energies and 
s.p.l. densities which determine the PSDs including the continuum effect [Figs. 9(b) and 
9(d)], with the related quantities for the particle-hole bound states [Figs. 9(a) and 9(c)]. 
The reasons for deviations from the PSD general trend, for the two main additional classes 
of particle-hole state densities with variant characteristics are outlined in this way. Thus, 
the average excitation energies, limited by the value B for the excited particles, entirely 
determine the s.p.l. densities for the particle-hole bound state densities. Actually one may 
note that g p has a rather constant value in this case. The value of gn is similarly constant but 
obviously lower in the latter case, i.e. of the bound and continuum state density including 
the continuum effect. However, the above-discussed aspect of g p is the key quantity for the 
second class of modified PSDs. 

Moreover, these two classes of PSDs with various restrictions [Figs. 10(a) and 10(b)] 
are at the same time compared with the predictions of Eqs. ([23|) and (|24[) including only 
the finite-depth correction, also shown in Fig. 10(b). A few configurations significant in 
PE calculations are used in this respect. The ratios of each of the two variant PSDs to 
the general PSD values given in Fig. 10(b) are further shown in Figs. 10(c) and 10(d). It 
follows that at medium energies the size of the continuum effect on the PSD values, which 
is given by the latter class of ratios, is rather similar to that for the bound-state condition. 
Therefore, a possible replacement within MSC calculations of the particle-hole bound state 



density by the PSD corrected for the continuum effect [27J] would not be trivial. A similar 
point may concern the use within the MSD calculations of the PSD either including the 
continuum effect PSJ or taking into account the free-gas single particle levels as well. 



The relation between the results of Eq. (^) and the quantum-mechanical (QM) calcu- 
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lations should be also considered before further use of the former in reaction calculations. 
Hence, following Shlomo et al. fl5C| , |SflJ| , the s.p.l. density was calculated by using the 



respective relation with Green's function. As an alternative to the smearing procedure, the 
imaginary part of Green's function has been calculated separately for the discrete and con- 
tinuous states. The regular and Jost solutions of the radial Schrodinger equation are used 
in the continuum. The smooth part of the rapidly fluctuating s.p.l. density is calculated by 
means of the Strutinski smoothing procedure j49|. The Woods-Saxon (WS) potential |30fl 
was considered in this frame as well as in the semiclassical Thomas-Fermi (TF) formula, 
with the similar results shown in Fig. 11(a). The familiar FGM shape is given by the 
TF formula with an infinite square-well (SQ) potential, while the corresponding finite well 
(FSQ) illustrates the continuum effect in Fig. 11(b). It should be noted that the continuum 
component of the s.p.l. density is nearly the same within either exact quantum-mechanical 
calculations with the WS potential, or TF approximation with either the WS or the FSQ 
potential wells, provided that the free-gas contribution is subtracted. Moreover, a similar 
trend is obtained by means of the simple FGM formula Q25D taking into account the con- 
tinuum effect. Nevertheless, the quantum-mechanical s.p.l. density can be related to this 
formula only for a reduced Fermi energy, e.g. F ~ 20 MeV [|50|| . This value has been ob- 



tained as an average value along the trajectory of the incident projectile with respect to the 
both nuclear density and first nucleon-nucleon collision probability. The usual value F=38 
MeV causes lower g(sp) values, that are not consistent with the phenomenological data. 

Therefore one may use the simple FGM energy dependence, within an appropriate form 
which matches the quantum-mechanical s.p.l. density including the continuum effect, in 
the average energy-dependent ESM formalism. This unsophisticated yet improved method 
could provide the correct PSDs for MSD/MSC calculations, in agreement with the consid- 
eration that the highly-excited single-particle states are not strongly coupled to compound 



nuclear states |oT| or partially relaxed states of composite nuclei formed in nuclear reactions 
at intermediate energies |52] , |53| . The question is additionally made intricate by the recent 
proof of a much shorter time scale required to reach thermal equilibration in intermediate- 



energy nucleon-induced reactions, found to be of the order of ~10 -22 sec ||54j| . Further 
experimental-data analyses should thus consider a combination of reaction models and re- 
lated PSD formalisms as well. 

V. SUMMARY AND CONCLUSIONS 

The particle-hole state density has been obtained by means of recursive relations, for the 
bound as well as bound-plus-continuum states. The corresponding expressions, i.e., Eqs. 
(|) and (p|), can be used for various energy dependences of the excited-particle and hole 
state densities while the particular case of the FGM is discussed. We have underlined that 
consideration of the finite depth of the nuclear-potential well should be explicitly present in 
the particle-hole state density formulas exceeding the ESM framework. 



Next, the results of the recursive formula are compared with the Kalbach |L8| approxi- 
mation still within the ESM formula but using distinct average s.p.l. densities for the holes 
and excited particles respectively, at their average excitation energies. At the same time 
the Kalbach formalism is extended to the case of the bound states, while the pairing and 
Pauli-blo eking effects have been included in all terms of the ESM correction function. The 
correctness of the average energy-dependent ESM approach is established by reference to 
the rigorous convolution (H) of the energy- dependent s.p.l. densities for the case of the 
FGM dependence. The difference between the predictions of the two methods is compared 
with the similar variation between average energy-dependent form and the standard formula 



35| , |36|1 , the former being much lower especially in the bound-state case. 

The exact calculation of the Pauli-blocking effect, which is close to the well-known ap- 
proximate Pauli correction |5|,[J!|[36|] , is also discussed. Thus it is shown on a common basis 



the higher importance of the s.p.l.-density energy dependence versus the exact calculation 
of the Pauli correction. The significant role of the pairing correction is pointed out, while 
comments are made on the circumstances under which the adequate account of the pairing 
effects could indeed appear less than essential HH . 
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The continuum effect has been considered for the case of a FGM energy dependence 
in the average energy- dependent ESM approach. The continuum component of the s.p.l. 
density is found rather similar using either exact quantum-mechanical calculations with the 
Woods-Saxon potential, or Thomas-Fermi approximation with WS as well as finite-square 
potential wells, provided that the free-gas contribution is subtracted. A similar trend is 
obtained by means of the simple FGM formula for the s.p.l. density if the continuum effect 
is taken into account. It should be noted that no arbitrary truncation, e.g. in the range 
15-25 MeV |38|| , is thus necessary in order to take care for the continuum effect within the 
s.p.l. density account. On the other hand, since the actual quantum-statistical analyses 
of the multistep reactions use the rough ESM, the results following consideration of the 
effective AA-interaction strength as the only free parameter could be altered. This point 
is the subject of current work along with systematic calculations of the s.p.l. density in the 



continuum and the correlation with PE surface effects |18|,|50|. 
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APPENDIX: FERMI GAS MODEL S.P.L. DENSITY IN CONTINUUM 

The s.p.l. density associated with a local mean field V has the following expression in the 
Thomas-Fermi approximation, by taking into account the spin degeneracy and neglecting 
the spin-orbit interactions p8|.|30|.[55| 
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9 TF (e) 



1 /2m\ 3/2 



J dr (e- V(r)) 1/2 Q(e- V(t)) . 



(Al) 



The single-particle energy e is measured relative to the top of the nuclear well, in order to 
make a clear distinction between the bound states at e <0 and unbound continuum states 
at e >0. 

For the finite well potentials, the nucleus can be imagined inside a spherical box of radius 
R larger than the range of V(r) (see Fig. 1 of |HJ). In the case of a square potential well of 
radius R and depth V <0, we have from Eq. (JAT1) 



9l F {e) 



/2m\ 3/2 r 



n (e - y ) 1/2 e( e - v ) + fie^efc) - cV /2 ©(e)] , (A2) 



2tt 2 V h 2 J 

where f2 =47ri?g/3 and Q=4iiR 3 /3. Since the properties of the nucleus itself do not depend 
on R 1 2l| , in the case of the finite square well one has to subtract the contribution of free 
Fermi gas when e > [ ^ , PH| ,PT[] 

Q /2m\ 3/2 



,1/2 



(A3) 



The s.p.l. density which is thus obtained 

1 /2m\ 3 / 2 



9fsq( ( 



2tt 2 V Ti 2 J 



v )^e{e-v )-^e(e)] 



(A4) 



has the well-known FGM form except the continuum correction term. In terms of the 
single-particle energy e=e-Vo which is measured relative to the bottom of the nuclear well, 
it becomes 

e\ 1 / 2 re + V \ 1/2 



9fsq(* 



9o 



Q(e + V ) 



(A5) 



where go=g(eF)=g(F) with reference to both notations used for the s.p.l. energy. Actually, 
the derivation of Eq. (|A4|) shows that the radius R is indeed taken into account but gpsoi 6 ) 
does not depend on it. Therefore, the final expression is apparently only a difference of 
terms calculated for the SQ potential well and the free-particle case, respectively, within an 
infinite spherical box with the radius Rq. 
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FIGURE CAPTIONS 



FIG. 1. The convolution state density given by Eq. (Q), divided by the Oblozinsky formula 
without Pauli correction (solid curves), as well as by the Ericson formula (dashed 
curves), for the given p-h configurations. In Eq. (|9]) and the Oblozinsky formula, 
F=38 MeV and the limit of large B is considered. For all calculations go=14 MeV -1 . 

FIG. 2. The correction function for the Pauli blocking, nuclear-potential finite depth and 
pairing-correlation effects, of the ESM particle-hole state density for given p-h con- 
figurations. The calculations use the finite-depth correction alone with F=38 MeV 
(dotted curves), the Pauli correction without (dashed curves) and with the pairing 
effects included (solid curves), and also the bound-state condition with 5=10 MeV 
(dash-dotted curves). 

FIG. 3. (a),(b) The f^-(p,h,E,F) and (c),(d) the f K (p,h,E,F) correction functions to the 
ESM formula for the Pauli blocking, potential finite-depth and pairing-correlation 
effects, as well as (e),(f) the average excitation energies for particles and holes within 
the given p-h configurations, for (a),(c),(e) the particle- hole bound plus continuum 
states, and (b),(d),(f) the bound states only. For all calculations F=38 MeV, and 
B=10 MeV is considered for the bound states. 

FIG. 4. (a) The h, E, F) and (b) Jk{p, h, E, F) correction functions to the ESM formula 
for the Pauli blocking, potential finite-depth and pairing-correlation effects, (c),(d) 
the average excitation energies for particles and holes, and (e),(f) the average energy- 
dependent s.p.l. densities for excited particles and holes, for the 2plh configuration. 
The meaning of the curves is the same as in Fig. 2 except that the dotted line in (e) 
and (f) gives the value of go. For all calculations go=14 MeV -1 and F=38 MeV, while 
B=10 MeV is considered for the bound states. 

FIG. 5. The particle-hole state densities for the given p-h configurations, obtained with the 
average energy-dependent ESM formalism (solid curves) and the PSD recursive formula 
(dashed curves), and their ratio, for (a),(c) bound-plus-continuum states, and (b),(d) 
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bound states only. The same global values are used as in Fig. 4. 

FIG. 6. The same as Fig. 5 except the latter PSD formula considered in the comparison is the 
ESM formula (dash-dotted curves). 

FIG. 7. (a),(b) The particle-hole state densities for given n-exciton configurations with p—h, 
obtained with the ESM formulas including the advanced pairing correction jTT2| and 
either the exact correction for the Pauli exclusion principle (solid curves) or the 
respective approximate form (dashed curves), and (c),(d) their ratios. For all 
calculations go=14 MeV -1 , while Ao=l MeV is considered for the pairing correction 
account in (c) and (d). 

FIG. 8. Average excitation energies of the excited particles (dotted curves) and holes (dashed 
curves) for given n-exciton configurations with p=h, as function of the total excitation 
energy. The same global values are used as in Fig. 4, while A =l MeV is considered 
for the pairing correction account. 
FIG. 9. (a),(b) The average excitation energies and (c),(d) the related s.p.l. densities for (a),(c) 
the particle- hole bound states, and (b),(d) the general case including (solid curves) or 
not (dotted curve) the continuum effect, for the 2plh configuration. The same global 
values are used as in Fig. 4. 

FIG. 10. The particle-hole state densities for the given p-h configurations, obtained within the 
average energy-dependent ESM formalism for (a) the bound states and (b) the bound- 
plus-continuum states including the continuum effect (solid curves) or taking into 
account also the free-gas single particle levels as given by Eqs. fl23|) and (|24"D in the 
limit of large B and with only the finite-depth correction (dotted curves), and (c),(d) 
the ratios of each of the two kinds of variant PSDs to the third one. The same global 
values are used as in Fig. 4. 

FIG. 11. The comparison of the smoothed quantum-mechanical s.p.l. density for the neutrons 
of the nucleus 56 Fe (solid curve) and the results of the TF approximation using (a) the 
same Woods-Saxon potential well as within the QM calculation, and (b) the infinite 
(SQ) and finite square potential wells (FSQ) (dashed curves), and the FGM formula 
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with the Fermi-energy values of F=38 MeV (dotted curve) and F=20 MeV (dot-dashed 
curve). For parameters of the potential wells see Refs. [pOj , ^] . 
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